INSTABILITY OF STEADY STATES FOR NONLINEAR WAVE AND HEAT 

EQUATIONS 



PASCHALIS KARAGEORGIS AND WALTER A. STRAUSS 

Abstract. We consider time-independent solutions of hyperbolic equations such as duu — 
Alt = f{x,u) where / is convex in u. We prove that linear instability with a positive eigen- 
function implies nonlinear instability. In some cases the instability occurs as a blow up in finite 
time. We prove the same result for parabolic equations such as dtu — Au = f{x, u). Then we 
treat several examples under very sharp conditions, including equations with potential terms 
and equations with supercritical nonlinearities. 



1. Introduction 

Given a linear second-order elliptic differential operator L whose coefficients are smooth and 
bounded, consider the parabolic equation 

dtu + Lu = f{x,u), xgM" (1.1) 

along with its hyperbolic analogue 

d'fu + adtU + Lu = f{x, u), x e M", (1.2) 

where / is a nonlinear term and a G M is arbitrary (possibly zero). A very important step in 
understanding the behavior of general solutions lies in understanding the qualitative properties 
of special types of solutions. In this paper, we focus on time-independent solutions, also known 
as steady states, and we address their stability properties in the context of both (11. ip and (11. 2p . 
Our main goal is to provide sufficient conditions under which linearized instability can be used 
to draw conclusions about nonlinear instability. 

Before we turn to our main results, however, let us first introduce some assumptions on the 
time-independent solution ip and the nonlinear term /. We are going to assume that 

(Al) the equation Lip = f{x, ip) has a solution ip; 

(A2) the adjoint linearized operator L* — fu{x,ip) has a negative eigenvalue — and a cor- 
responding eigenfunction x ^ /^^(M") fl L^(M") that is non-negative; 
(A3) both f{x,ip) and fu{x,ip) are bounded; 
(A4) the nonlinear term f{x,s) is convex in s and is C^. 

Here, (A2) is mostly meant to ensure the presence of a negative eigenvalue; our assertions about 
the eigenfunction x are already implied by this fact under pretty general conditions. Also, note 
that we do not require the steady state ip to be bounded. 

A non-technical description of our main result is that if is a nonlinearly unstable solution 
of both (II. ip and (11.21) whenever (Al)-(A4) hold; namely, solutions which start out close to ip 
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need not remain close to ip for all times. When it comes to the parabolic case (11.11) . our result 
applies to a wide class of initial data, including all initial data for which 

u{x,{])> ^{x). (1.3) 

To establish the instability of we show that the norm 

||u(a;,t) - v2(x)||l°c(R") (1.4) 

must either grow exponentially at all times or else blow up in finite time. Our main result 
for the hyperbolic case (11. 2p is almost identical. It too applies to a wide class of initial data, 
including all initial data for which 

m(x, 0) > (^(x), dtu{x, 0) > (1.5) 

and it shows that the energy norm0 

\\u - ip\\e = \\u{x,t) - ip{x)\\HHR^) + \\dtU{x,t)\\L-i{K.^) (1.6) 

must either grow exponentially at all times or else blow up in finite time. 

Note that the above results prove instability in a sense that is much stronger, and perhaps 
more natural, than the usual one. Namely, not only do they show that the solution exits any 
given neighborhood of in finite time, but they also ensure that the solution does not reenter 
that neighborhood at any later time. 

If one is willing to impose an additional positivity condition on the nonlinear term, then the 
above results can be further improved to show that instability occurs by blow up. This is the 
case, for instance, if one additionally assumes that 

(A5) there exist Co > and p> I such that /(x, s) > Cq\s\p for all (x, s) G M" x M; 
(A6) the product ipx is integrable, where x denotes the eigenfunction from (A2). 

In our main results, we shall prove instability assuming (Al)-(A4), and also instability by blow 
up assuming (Al)-(A6). As the reader should note, (A3)-(A6) may be safely ignored in the 
special case f{u) = \u\^ for some p > 1, provided that ip is bounded. For that special case, in 
particular, we establish instability by blow up assuming (Al)-(A2) only. 

To a large extent, our results are complementary to the abstract result of Grillakis, Shatah 
and Strauss [9j. These authors deal with arbitrary nonlinearities and the more general class of 
bound states, but their assumptions are more restrictive than ours. As we already mentioned, 
another advantage of our approach is that our notion of instability is much stronger than the 
one used in [9]. As we indicate later in the introduction, a final advantage is that our result 
applies to a very sharp class of perturbations, which is not the case with the result of [9]. 

Section [2] is devoted to the proof of our two main instability results. Theorems 12.11 and 12.21 
Our proof is quite elementary and based on a variant of Kaplan's eigenfunction method [TT] 
which allows us to reduce our analysis to the study of a certain functional. Although there are 
several other variants of this method, they are concerned with the associated Cauchy problem 
rather than the instability problem; see [3, [8], [121 [HI [13 [27], for instance. In the last three 
sections of this paper, we apply our two main results to treat some examples that we briefly 
discuss below. 



In each case, our norm agrees with the one dictated by the local existence theory. 
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Our first and most general example appears in Theorem 13. 1[ Here, we shall only describe a 
very special case of the theorem that is itself of independent interest. Consider the nonlinear 
heat and wave equations with potential! 

diu- Au + V{x) ■u = \u\P, xeW; i = l,2. (1.7) 

We take p > 1 and assume that V is bounded, continuous and non-negative, but we make no 
other assumptions. Then Theorem 13.11 implies the instability of all non-negative, steady 
states that vanish at infinity. A sufficient condition on V which ensures the existence of such 
steady states (for some p) is provided by [22] and stated in Theorem l3.3[ Together with Theorem 
13. H these results impose no restrictions on the linearized operator; we are not aware of any 
other general results with this feature. 

In our second example. Theorem 14.21 we focus on the two-dimensional equation^ 

diu - Am = e", xE M^; i = l,2. (1.8) 

The classification of all steady states for which the nonlinear term is integrable is provided 
by [1] and stated in Lemma 14.11 They are unbounded and, as far as we know, their instability 
was not previously known. Theorem 14.21 shows they are all unstable. 
In our last example. Theorem 15. 6^ we focus on the equations 

diu- Au=\u\P, xeW; i = l,2, (1.9) 

where n > 2 and p > When it comes to the existence of positive, C^, radially symmetric 
steady states, the known results [31 HI El [IT] are summarized in Theorem 15. 2[ Except for the 
special case p = these steady states are not in H^, so our first example is no longer 
applicable. In section [5l we analyze the spectrum of the linearized operator and we find the 
exact values of p for which a negative eigenvalue emerges. Using our main results, we then 
show that there is a critical value pc, depending on the space dimension, such that the steady 
states are nonlinearly unstable if p < pc and linearly stable if p > Pc] see Theorem 15.61 

Our conclusions for equation fll.9p are new only in the hyperbolic case. As for the parabolic 
equation, the critical value Pc emerged in [TUj, where a quite different approach was used. The 
instability result of [lOj is weaker than ours because it applies to a smaller class of initial data, 
but the stability result given there is stronger since it proves nonlinear stability instead. On 
the other hand, the approach in [10] relies on the maximum principle, so it cannot be applied 
to yield analogous results for the wave equation; see Remark 15.71 for more comments. 

Finally, our result for the wave equation (11.91) is closely related to a recent paper of Krieger 
and Schlag [13]. These authors focus on the three-dimensional quintic case 

d^u- Au=\u\\ xeR^ (1.10) 

and a particular positive, steady state (p. Theorem 15.61 implies that ip is unstable, while the 
result of [13] asserts the existence of a stable manifold associated with ip. Based on numerical 
evidence from [21I21I, this stable manifold ought to separate the small perturbations for which 
solutions exist globally from those for which solutions blow up. In fact. Theorem 15.61 provides a 
partial proof of this conjecture, namely the blow up for all perturbations that lie strictly above 
the tangent plane to the stable manifold at the origin; see Remark 15.81 for more details. 



^In the rest of the introduction, the possible damping term adfU is suppressed merely for the sake of exposition. 
''This is our only example for which the space dimension is restricted. 
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2. The main instability results 

In this section, we prove our two main instability results regarding the solutions of a general 
elliptic equation of the form 

L<^ = f{x,<^), xeW. 

Recall that L is a linear, second-order elliptic differential operator whose coefficients are smooth 
and bounded, while / is a nonlinear term. Our precise assumptions (Al)-(A6) were already 
mentioned in the introduction, so we shall not bother to repeat them here. 
First, we deal with the hyperbolic case and thus focus on the equation 

dfU + adfU + Lu = f{x,u), u{x,0) = {p{x) + ipo{x), dtu{x,0) = ipi{x). (2.1) 

Since the steady state ip is an exact solution when ipo = ipi = 0, we are mostly concerned with 
the case that the perturbation (^/'q? V'l) is small in some sense. In our next result, we shall deal 
with finite-energy perturbations. Note, however, that ip itself need not be of finite-energy. 

Theorem 2.1. (Hyperbolic Equation) Let a G M. Assume (A1)-(A4) and let (^o,ipi) G 
i7i(M") X L2(M") be such that 



x{x) ipo{x) dx + / x(a;) ^/'i(x) > 0. (2.2) 
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Let < T < oo and let u be a solution of (12. ip on [0, T) such that u — (f is continuous in t 
with values in the energy space and f{x,u) is locally integrable. 

(a) IfT = oo, then the energy norm 

\\u{t) - IfWe = \\u{-,t) - V?(-)||h1(IR") + PtM(-,t)||L2(R") (2.3) 

must grow exponentially. 

(b) Assume also (A5)-(A6). ThenT < oo. 

Proof. Consider the function 

^i^) = / xi^) ' ^i^^^) dx, w{x,t) = u{x,t) — (p{x). (2.4) 



By (A2), one certainly has 

\G{t)\ < \\x\\lhr«) ■ |k(t)||L2(iR„) < C\\w{t)\U. 

Thus G(t) is well-defined and bounded as long as the energy remains bounded; the energy grows 
exponentially, provided that G{t) does; and the energy becomes infinite whenever G(t) does. 
In view of these facts, our assertions about the energy (12. 3p will follow once we have established 
analogous assertions for G{t). 

Let us first focus on part (a) and assume that w = u — ip is continuous on [0, oo) with values 
in the energy space. Then w is a solution of 

d^w + adfW + Lw = f{x, w + ip) — f{x, ip) 

in the sense of distributions. In view of our convexity assumption (A4), we have 

d^w + adtw + [L — fu{x, ip)] w >0. 
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Since > 0, we may multiply the inequality by x{^)(^{t)j where 6{t) is an arbitrary non- 
negative test function, and integrate to obtain 

xix) ■ dtw ■ 9'{t) dxdr + / ■ adtw ■ 9{t) dx dr 



+ [ [ [L* -Ux,ip)]x{x)-w-e{T)dxdT>0. (2.5) 
Jo Jr" 

To simplify the first two integrals, we note that 

G'{t) = / xix) ■ dtw{x,t) dx 



is a continuous function of t by our definition ( 12. 4p because dtW is continuous with values in 
L^(M"). To simplify the third integral in (12. 5p . we note that 

[L* - f^{x,ip)]x = -(T^X 
by our assumption (A2). Thus, equation (12.51) reduces to 

- [ G'{t) ■ 0'{t) dr + a [ G'{t) ■ ^(r) dr - / G{t) ■ 0{t) dr > 



Jo Jo Jo 

for all non-negative test functions 9{t), or equivalently, 

G"{t) + aG'it) - a^G{t) > (2.6) 

in the sense of distributions. Next, we note that our assumption (12. 2p reads 

^±^^^±^.G{0) + G'iO)>0. 

According to Lemma [2.31 below, both G{t) and G'{t) must then grow exponentially fast, so the 
proof of part (a) is complete. 

Next, we turn to part (b). Suppose that T = oo. As we have just shown, both Git) and 
G'{t) must grow exponentially fast. Using our assumptions (A3) and (A5), we have 

+ adtW + [L- /„(x, v?)] w = f{x, w + ip) - f{x, Lp) - /„(x, ip)w 

> Go\w + iff - Gi - Gi\w\. 

Multiplying by the non-negative eigenfunction x and integrating over space, we then easily find 
that 

G"it) + aG'it) - a^Git) > Go [ x\w + ^pf dx - Gi [ xdx-Gi [ xH dx, (2.7) 

the second derivative being understood in the sense of distributions. This can be justified, as 
above, by introducing a non- negative test function 6 it) to avoid the second-order derivative. 
Now, the last integral is at most 



by the triangle inequality. Since x ^ -^^ by (A2) and xf ^ by (A6), we deduce that 
G"it) + aG'it) > a^Git) + Go [ xk + dx - Gi [ x\w + (p\ dx - Cg. 
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Since G{t) grows exponentially fast, this implies 

G"{t) + aG'{t) >Co [ x\w + <^\^ dx-Ci [ x\w + </^| dx 

= CoA{t) - CiB{t) (2.8) 

in the sense of distributions for all large enough t. 

Now B{t) itself grows exponentially fast, as the triangle inequality gives 

Git)< I x\w\dx< I x\w + v\dx+ I xW\dx = B{t) + Cz. (2.9) 

JR" JR" JR" 

Similarly, A{t) grows exponentially fast since Holder's inequality gives 

B{t)<([ xdx] ([ xlw + ifl^dx] =C4A{ty/P. (2.10) 

This actually forces A(t) to grow faster than B{t), namely 

A{t)>C-^B{ty, (2.11) 

hence A{t) will eventually dominate B(t). Combining ([22]), (12^9]) and fl2TTTD . we now get 

G"{t) + aG'{t) > C^Ait) > G^Bitf > G^G{tY (2.12) 

in the sense of distributions for all large enough t. Moreover, both G{t) and G'{t) are eventually 
positive by above. Invoking Lemma [2.41 below, we reach the contradiction T < oo. ■ 

Next, we modify our previous approach to obtain a simple parabolic analogue of Theorem 

m 

Theorem 2.2. (Parabolic Equation) Assume (Al)-(A4) and letipo G L°°(M") be continuous 
with 



x{x)iJo{x) dx > 0. (2.13) 

Let < T < oo and let u be a solution of 

dtu + Lu = f{x,u), u{x,0) = ip{x) + ipo^x) (2-14) 
on [0, T) such that u — ip is continuous and bounded. 

(a) IfT = oo, then the norm \ \u — 9?||l°°(R") must grow exponentially. 

(b) Assume also (A5)-(A6). Then T < oo. 

Proof. Once again, it suffices to establish analogous assertions for the function 
^(^) = / xix) ■ w{x,t) dx, w{x,t) = u{x,t) — ip{x). 



Since 

\G{t)\ < ||x||l1{R") ■ ||'W||Lcx,(]Kn), 

the L°°-norm has to either grow exponentially or blow up whenever G{t) does. 

We now apply the same argument with minor changes. Note that w = u — ip satisfies 

dtw + Lw = f{x, w + ip) - f{x, ip), w{x, 0) = tpo{x) 
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in the sense of distributions. Arguing as before, we then get the estimate 

G'{t) - a^Git) > (2.15) 

instead of (12.61) . Since G'(O) > by assumption, the exponential growth of G{t) follows directly. 
Under the additional assumptions (A5)-(A6), our previous approach yields the estimate 

G'{t)-a^G{t)>Co [ x\w + ^ydx-Cif xdx-CifxHdx 

instead of (12. 7p . Then the remaining part of our proof applies verbatim to give 

G'it) > GjG{tY 

instead of (I2.12p . in the sense of distributions for all large enough t. Since G{t) is positive, it 
is easy to deduce that T < oo, as needed. ■ 

Lemma 2.3. Let a G M and b > 0. Suppose y{t) is a function such that 

y" + ay'-by>0 

on some interval [0,T) in the sense of distributions. If 

^^i^fi5-y(0) + y'(0)>0. (2.16) 

then both y{t) and y'{t) must grow exponentially on [0,T). 

Proof. Let Ai < < A2 be the roots of the characteristic equation + aA — 6 = and set 
z = y' — Xiy. Then 

z' - X2Z = y" + ay' -by>0. 
Using the test function exp(— A2t) 6{t), where 6(t) is another test function, it follows that 



[ z{t) exp(-A2r) 0'(r) dr > 0. 
Jo 



Choosing 6{t) to be an approximation of the characteristic function of the interval (0,t), we 
easily deduce that z{t) > e^'^^z{0). Thus y' — Xiy > e^^^z{0), which implies that 

yit) > e^^'y{0) + ^-^(0). 

A2 — Ai 

Since Ai < < A2 by above, the exponential growth of y then follows, provided that 



;.(0) = y'(0) - Aiy(O) = y'(0) + ^— 1/(0) 

is positive. In view of our assumption (12.161) . the exponential growth of y thus follows. 
Next, we use a similar argument to find that w = y' — \2y satisfies 

w' — \iw = y" + ay' — by > 0, 

hence also w(t) > e^^*w{0). Since y(t) grows exponentially by above, the equation 

y'{t)>X2y{t) + e^M0) 
then forces y'{t) to grow exponentially fast as well because Ai < < A2. 
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Lemma 2.4. Let a G M, 6 > and p > 1. Suppose y{t) is a non-negative function such that 

y{T,) > 0, y'{T,) > 0, y"{t) + ay'{t) > hyitf 
on some interval [Ti,T2) in the sense of distributions. Then T2 < 00. 

Proof. Suppose first that y E C^. Then the case a = 1 is treated for instance in Proposition 
3.1 of J25j; the case a > is quite similar; and the case a < is much easier. If y is merely C^, 
one can simply repeat the same argument using test functions as in the preceding lemma; we 
omit the details. ■ 



3. Convex nonlinearity with potential term 

In this section, we apply our main results to study non-negative if ^ solutions of the equation 

- A^ + V{x)-^ = f{^), xgM". (3.1) 

Although our approach applies verbatim to the more general case f{x,(p), the assumptions of 
our next result would have to be suitably modified in that case; we shall only deal with (13.11) 
for the sake of simplicity. 

Theorem 3.1. Assume that the following conditions hold: 

(Bl) V is continuous and bounded on M"; 

(B2) the essential spectrum of —A -\-V is contained in [0, 00); 

(B3) equation (13. ip has a non-negative solution ip G that vanishes at infinity; 
(B4) / G C^(M) is convex with /(O) = /'(O) = and f{^p) is not identically zero. 

Then (f is an unstable solution of any of the equations 

dtu + (-A + V)u = f{u), d'^u + adtu + (-A + V)u = f{u), 

where a G M. More precisely, the conclusions of Theorems \ 2. 1\ ( a) and \2.2( a) remain valid. 

Remark 3.2. By Weyl's theorem, our spectral assumption (B2) essentially requires that V{x) 
not be negative as |x| — > 00. We do not know of any existence results regarding (13. ip for which 
this assumption is violated; there are several existence results [U |22l [23] in case it holds. 

Proof. We verify the assumptions (Al)-(A4) of Theorems 12.11 and 12.21 Note that the existence 
assumption (Al) holds by (B3), while the convexity assumption (A4) holds by (B4). Since (p is 
bounded by (B3), both f{ip) and /'(v?) are bounded as well, so the assumption (A3) also holds. 
To check the remaining assumption (A2), we consider the self-adjoint linearized operator 

^ = -^ + v-n^). 

Noting that is a solution of (13. ip . we use (B4) to find that 

= /(^) - ^ ■ /'(^) = /(^) - /(O) - ^ ■ /'(^) < 0. 
Let ( , ) denote the standard inner product on L^(R"). Then the expression 

{^ip,ip)= f \v^\^dx+f [v-n^)]-^^dx 
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is non-positive because S^^p < < and is finite because Lp G H"^ and V — f'{f) G Thus, 
the first eigenvalue of ^ must also be non-positive. It can only be zero if 

is identically zero, in which case f{y) = yf'{y) for all y in the image of ip. Solving this ordinary 
differential equation, we get f{y) = Cy for some constant C. Then the fact that /'(O) = 
implies f{y) = for all y in the image of ip. Since this violates our assumption (B4), however, 
the first eigenvalue must actually be negative. 

Now that the first eigenvalue is known to be negative, the remaining assertions of (A2) follow 
by standard facts. Using a variational argument, it is well-known that the first eigenfunction can 
be chosen to be positive (see section 11 in [18], for instance). Finally, the first eigenfunction is in 
flL^ because it decays exponentially by Agmon's estimate; see Theorem C.3.5 in [21] or else 
a more general result of Nakamura [T9j. In fact, Agmon's estimate applies to any eigenfunction 
for which the associated eigenvalue lies below the bottom of the essential spectrum. Since 

hm /'(^(x)) = /'(O) = 
by (B3) and (B4), an application of Weyl's theorem (see Theorem 1.1 in [H]) gives 

(Tess(^) = aess(-A + V) C [0, Oo) 

because of (B2). Thus, the first eigenfunction must decay exponentially, as needed. ■ 

In our next result, we give a typical application of Theorem 13.11 Our precise assumptions 
on the potential V{x) are taken from [22]; they are merely meant to ensure that equation (13.31) 
has a non-negative solution ip G that vanishes at infinity. Needless to say. Theorem 13.11 
implies the instability of such solutions for a much wider class of potentials. The only advantage 
of our next result is that all the assumptions are verified explicitly. 

Theorem 3.3. Let n > 2 and 1 < p < Assume that V{x) is radially symmetric, continu- 
ous and locally Holder continuous with 

Ciil + \x\)-' <Vix) <C2 (3.2) 

for some constants Ci, C2 > and some < I < '^^'^"^^^'^^ ■ Then the equation 

-Aip + V{x) ■ <^ = <^P, xeW (3.3) 

has a positive, C"^ , radially symmetric and exponentially decaying solution ip. Moreover, ip is 
an unstable solution of any of the equations 

dtu + (-A + V)u = \u\P, d^u + adtu + (-A + V)u = \u\p, 

where a G M. More precisely, the conclusions of Theorems \2.1\ and \2.2\ (both (a) and (b)) remain 
valid. That is, the instability occurs by blow up in finite time. 

Proof. Under the given hypotheses, our assertions about the elliptic equation (13. 3p follow from 
Theorem 1.2 in [22J. These ensure that the assumption (B3) from the previous theorem holds 
with f{u) = \u\P. It is clear that this nonlinear term satisfies (B4). Moreover, V{x) is bounded 
and non-negative by (13.21) . so the remaining assumptions (B1)-(B2) hold as well. According to 
the previous theorem then, the conclusions of Theorems 12.11 and 12.21 remain valid, indeed. 
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To show that instabihty occurs by blow up, it remains to check that the additional assump- 
tions (A5)-(A6) of these theorems remain valid as well. Since f{u) = \u\^ satisfies (A5) and (p 
is bounded, these two assumptions hold trivially and the proof is complete. I 

4. Exponential nonlinearity in 2D 

In this section, we apply our main results to study solutions of the equation 

-A<f = e^, e^eL\R^). (4.1) 

The classification of all solutions is due to Chen and Li |i4J. It is perhaps worth noting that 
these solutions are unbounded and not necessarily of one sign. 

Lemma 4.1. Every C"^ solution of (14. ip is of the form 

ip{x) = log 32A^ ■ (4 + \^\x - yf)-^ 
for some A > and some ?/ G M^. 

Theorem 4.2. Every solution of (14. ip is an unstable solution of any of the equations 

dtu — Am = e", d^u + adtu — Au = e", 
where a G M. More precisely, the conclusions of Theorems {Kj^ a) and \2.2( a) remain valid. 

Proof. We set f{u) = e" and verify the assumptions (A1)-(A4) of Theorems 12.11 and 12.21 
The existence assumption (Al) and the convexity assumption (A4) obviously hold in this case. 
Assumption (A3) holds as well because 

f{^) = f'i^)=32\'.iA + \'\x-y\'r' 

is bounded. To show that (A2) holds as well, we focus on the linearized operator 

^ = -A -expy;(x) = -A-32X'^- {4 + X'^\x - y\'^)~\ 

Its essential spectrum is 

(Tcss(=^) = crcss(-A) = [0, oo) 
by Weyl's theorem since f'{^p) is bounded and vanishes at infinity. Exactly as in the proof of 
Theorem 13.11 the first eigenfunction will be positive and exponentially decaying, so long as the 
first eigenvalue is negative. In particular, it suffices to show that the associated energy 

^(0=/ |VC(x)|2rfx- / exp^(x)-C(x)'rfx (4.2) 

is negative for some test function ( G H^(R'^). Choosing 

ax) = {4 + X'\x-y\y\ 

the energy (14.20 is given by 

dx 





X — 


y 


^ dx 


(4 + A2 


X 


-y 


2y 



E{C) = 16A^ / t^^tP ^ - 32A 



(4 + A2|x-i/|2)6 



32AV / —^-64X^71 



(4 + A2r2)6 (4 + A2r2)6 
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and then a short computation gives 



TT 



/»00 /»CXD 

E{Q = IQn / (s - 4)s-^ ds - 32tt / s'^ ds 

J A J 4 



320 

This imphes the existence of a negative eigenvalue and completes the proof. ■ 

5. Power nonlinearity with zero mass 
In this last section, we focus on positive solutions of the equation 

- Aifix) = ifiixY, xeW. (5.1) 
Before we state our results, it is convenient to introduce the quadratic polynomial 

Q{a) = a{n-2-a). (5.2) 
This polynomial arises naturally through the computation 

-A|x|"" = a{n - 2 - a) ■ \x\~''-^ 
and it is closely related to Hardy's inequality 

Vu{x)\'^ dx> (^^^—^^ [ \x\'^u{xf dx, (5.3) 

which is valid for all u G //^(M") and each n > 2. Namely, the coefficient on the right side of 
(15. 3p is the maximum value of Q{a) and it is known to be sharp in the following sense. 

Lemma 5.1. Let n > 2 and let V be a bounded function on which vanishes at infinity. If 

there exists some e > such that 

V{x)<-{l + e)-(^^J -Ixr 

for all large enough \x\, then the operator —A + V has infinite negative spectrum. 

For a proof of Hardy's inequality (15. 3p . see pg. 169 in [20]. For a proof of the lemma, see the 
appendix in [[5], for instance. 

We are now ready to state the known existence results for positive steady states. Parts (a) 
and (b) can be found in either |3] or |3]. For part (c), see Theorem 1 in [17] and Theorem 5.26 
in [6]. 

Theorem 5.2. Let n>l and p > 1. Denote by Q the quadratic in (15. 2p . 

(a) // either n = 1,2 or p < then equation ( 15. ip has no positive, solutions. 

(b) If n > 2 and p = then any positive solution of (15. ip is of the form 



+ |x — y\ 



for some A > and some y G M". 
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(c) If n > 2 and p > then the positive, C"^ , radially symmetric solutions of (15.11) form 
a one-parameter family {v5a}a>0j where each (fa satisfies 

(^,(0) = a, lim |xp ^aixf-^ = Q f > 0. (5.5) 

|x|^oo \p — I J 

In what follows, we use the previous two results to analyze the spectrum of the linearized 
operator associated with (15. ip . In particular, we find the exact values of p for which a negative 
eigenvalue emerges. Although a preliminary characterization of these values is provided by the 
next two lemmas, a more concrete characterization will be given in Theorem 15. 6[ 



Lemma 5.3. Let n > 2 and p > Let Q be the quadratic in (15. 2p and let (p denote any one 
of the steady states provided by Theorem \5.B. Then —A — pip'^'^ has a negative eigenvalue if 

In particular, it has a negative eigenvalue if p = ^^f. 

Proof. Suppose first that p > Then (15. 5p and (15.60 allow us to find some e > such that 



lim 1x1' ^{xY = Q >p^\l + 2e) 

\x\^oo \p — 1 



n — 2 



Thus 

V{x) = -p<p{xY-^ < -(1 + e) ■ [^-^^j ■ 1^ 

for all large enough so the existence of a negative eigenvalue follows by Lemma [5. 1[ 
Suppose now that p = Then inequality (15.60 automatically holds because 



p — 1 J \p — 1 

for this particular case. According to part (b) of Theorem 15.21 we also have 

-A-pip{xy-^ = -A-X^n{n + 2) ■ [X^ + \x - y\^y^ 
for some A > and some y G M". Thus it suffices to check that the associated energy 

^(0=/ \VCix)\Ux- [ pv{xr-\{xfdx (5.7) 



is negative for some test function ^ G H^{W^). Let us now consider the test function 

C(a:) = (A2 + |x-yn-"/^-\ 

whose energy (15. 7p is given by 



dx 



+4 



„ {X^ + \x - y\^)^- 

Except for a positive factor, this expression is equal to 

~, , , , r°° r"+i dr r r-^-^dr , , 

m-in + 2)f ^-,_^_A^,| Jy^^.- (5-8) 
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Moreover, an integration by parts gives 

roc ^n+l ^ roc 



(A2 + r2)"+4 2{n + 3)Jo (A2 + r2)"+3 



n 



°° r""-^ dr ^ r"+i dr 



(A2 + r2)"+4 7o (A2 + r 



2{n + 3) 
hence also 



2^n+4 



(A2 + r2)«+4 n + 6Jo (A^ + r2)"+4 
Inserting this equahty in (15.81) . we arrive at 



+ 6 y Jo (A2 + r2)"+4 n + 6 (A2 + r2)«+4 ' 
Since this expression is negative, the energy (15.71) must also be negative, as needed. 

Lemma 5.4. Let n> 2 and p > Let Q denote the quadratic in (15. 2p and assume that 



n-2V „ / 2 



>P-Q\ 7 • (5.9) 



2 J \p-l 
Then the steady states provided by Theorem \5.2\ satisfy 

\x\^ ^{xy-' < Q (^-^^ (5.10) 

for each x G M", and the operator —A — pip'^^^ has no negative spectrum. 

Proof. First suppose the estimate (I5.10p does hold. Using our assumption (15.91) . we then get 

-p^ixf-' >-p-Q ■ \xr > - (^) 

for each x G M". According to Hardy's inequality (15.31) . this already implies that 

\Vu{x)\'^ dx — / p(p{xY~'^ ■ u{xy dx > 



2 

■ \x\ 



for all u G H^{W^), so that the operator —A — pip'^ ^ has no negative spectrum, indeed. 
Let us now prove ( ]5.10p . Set k = and consider the function 

W{s) = e^^'ipi^e^) = r^ip{r), s = logr = log 
Then W{s) is positive and bounded because of (15. 5p . It also satisfies the equation 

Q{k-ds)W{s) = W{sY (5.11) 
by Lemma [5.51 below. We note that s ranges over (— oo, oo) as r ranges from to oo, while 

lim W{s) = lim r^p{r) = 0. 

The derivatives of W{s) must also vanish at s = — oo because 

lim W'{s) = lim r ■ dr[r'^ip{r)] = 0, 

s— >— oo f^o+ 
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and so on. Now the fact that f{x) = is convex on (0, oo) ensures that 

fiwis)) - f (g(fc)^) > /' (Q(fc)^) ■ (wis) - g(A:)^ 

that is, 

W{sY - Q{k)^^ > pQ{k) ■ (W{s) - Q(fc)^) . 
Inserting this inequahty into ( 15. lip , we thus find 

Q{k - d,)Wis) - Qik)^ > pQ{k) ■ {w{s) - Qik)^^ . (5.12) 
To ehminate the constant term on the left hand side, we change variables by 

Y{s) = W{s) -Q{k)^. 
Then equation (15.121) can also be written in the equivalent form 

pQik)-Q{k-d,)^Yis) <0. (5.13) 

We will prove below that the characteristic polynomial 

^{X)=pQ{k)-Q{k-X) (5.14) 

has two negative roots Ai, A2. Temporarily assuming this and noting that ^{X) has its highest- 
order coefficient equal to 1 by (15.21) . we can then factor the ordinary differential inequality 
(I57[3|) as 

Z{s) = Y'{s) - XiY{s), Z\s) - A2^(s) < 0. 
The last inequality makes e~^'^'^Z{^s) decreasing, so we must actually have 

e-^^'Z{s) < lim e~^''Z{s) = 

because A2 < by above. This gives Y'{s) — XiY{s) = Z{s) < 0, and we similarly find 

e-^''Y{s) < lim e-^''Y{s) = 

because Ai < as well. Now that Y{s) < 0, we have 

W{s) < Q{k)^. (5.15) 

2 

This is precisely the desired inequality (15.101) since W{s) = rp-^ip{r) by definition. 

Thus it remains to show that the quadratic in (I5.14p has two negative roots in case of 
assumption (15. 9p . Since the product of the two roots is 

^(0) = (j9 - 1) • g(A;) > (5.16) 

by (I5.14p and (15. 5p . it suffices to show that ^{X^) < for some A,,, < 0. We choose 

, n-2 2 n-2 

X^ = k = 

2 p-l 2 

and note that A* < because p > by assumption. Combining (I5.14p . (15. 2p and (15. 9p . we 
find ^ 

j^{X,)=pQ{k)-Q(^)=pQ{k)-(^) <0. 
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Thus ^{X) has two negative roots. ■ 

Lemma 5.5. Let n > 2. Let p > and k = Given any one of the steady states if 

provided by Theorem \5.S\ the function 

W{s) = e''V(e'') = r''{p{r), s = \ogr = log 

must then satisfy the ordinary differential equation 

{ds -k){ds-k + n-2) W{s) = -W{sy. (5.17) 

Moreover, this ordinary differential equation can be written in the form 

Q{k-ds)W{s) = W{sY, 

where Q denotes the quadratic in ( \5.2\\ . 

Proof. Since dr = e~^ds, a short computation allows us to write the radial Laplacian as 

A = d^ + {n- l)r-^dr = e-^'{n -2 + ds)ds. 
Since —Aip = ip^ and kp = k + 2hy assumption, this easily leads us to 

-lyp = -e^P'^P = e^^+2^A(^ = e'''{n-2 + ds)ds{e-^'W). 
Using the operator identity dse~''^ = e~^^{ds — k) twice, we then get 

-WP = e'''{n -2 + ds)e-'"'{ds - k)W = {n-2 + ds- k){ds - k)W. 
This proves our first assertion fl5.17p . from which our second assertion follows trivially. ■ 

Theorem 5.6. Let n > 2 and p > Let ip denote any one of the steady states provided by 
Theorem \5.2\ and set 

( oo ^f n<10 \ 

Pc=< 8n+4+8Vp ^/ ^ > IQ f ^^'^^^ 

I (n— 2)(n— 10) ■' I 

for convenience. Then pc > and the following dichotomy holds. 

If p < Pc, then if is a nonlinearly unstable solution of any of the equations 

dtu -Au= \u\P, d^u + adtu - An = \uf, (5.19) 

where a G M. More precisely, the conclusions of Theorems \2.1\ and \2.2\ (both (a) and (b)) remain 
valid. That is, the instability occurs by blow up in finite time. 

Ifp > p^, on the other hand, then (p is a linearly stable solution in the sense that the linearized 
operator has no negative spectrum. 

Proof. Let Q denote the quadratic in (15.21) and consider the expression 



Q{p) =A{p-l 



\2 



n-2V „ / 2 



-P-Q 



p — 1 



(5.20) 



When this is non-negative. Lemma 15.41 implies that ip is linearly stable. When it is negative, 
on the other hand. Lemma 15.31 implies that the linearized operator has a negative eigenvalue. 
Since (p vanishes at infinity, each of the functions 
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is bounded, whence (A3)-(A6) hold and the linearized operator has essential spectrum 

aess(-A - pipP-^) = aess(-A) = [0, oo), 

whence (Al)-(A2) follow by standard arguments, as in the proof of Theorem 13. 1[ 
In view of these observations, it thus suffices to check that 

Q{p)<0 ^ !i±|<p<p^. (5.21) 

Combining our definitions (15.21) and (I5.20p . we have 

Q(p) = (n - 2)(n - 10) • p^ - 2{n^ - 8n + 4) ■ p + (n - 2)1 (5.22) 
The last equation easily leads to 

\n-2j n-2 ^ ^ 

Case 1. Suppose first that n > 10. Then Q{p) has two roots, the largest one of which is Pc, 
given by (I5.18p . Since Q{p) is negative only between its two roots, (I5.23P implies that Q{p) is 
negative if < p < Pc and non- negative if p > Pc- 
Case 2. Suppose now that 2 < n < 10. Then the computation 

Q'{p) = 2{n - 2)(n - 10) (^p - ^^^^ - 48 

forces Q{p) to be decreasing for each p > and hence negative by (I5.23p . This shows that 
the desired condition (I5.2ip holds with pc = oo. ■ 

Remark 5.7. When it comes to the nonlinear heat equation (I5.19p . the above result of Theorem 
I5.6l is originally due to Gui, Ni and Wang JTU\. Their proof is based on an earlier result of Wang 
[26j, according to which the graphs of two steady states do not intersect if p > Pc, while the 
graphs of any two steady states do intersect if p < pc- Using this characterization of the 
critical value Pc, they are then able to employ comparison arguments which rely on the strong 
maximum principle. As far as the unstable case p < Pc is concerned, they show that blow up 
occurs provided the perturbation ipo in (12.140 satisfies < ipo ^ 0, while solutions exist globally 
provided > V'o ^ 0. Theorem 15.61 yields a refinement of their blow up result for the nonlinear 
heat equation since our exact assumption (12.130 on ipo reads 

x{x) ipo{x) dx > 

for a certain positive function x{^)- As for the stable case p > pc, the result in [10] is much 
stronger than Theorem 15.61 because it proves nonlinear stability, not merely linear stability. 
However, a proof of nonlinear stability for the analogous hyperbolic problem is difficult to 
imagine. 

Remark 5.8. When it comes to the wave equation (I5.19P in the critical case p = ^^1, our con- 
clusions are closely related to a result of Krieger and Schlag [13] in three dimensions. According 
to Theorem 15.21 the equation 

d'^u -Au= \uf, X eW^ (5.24) 
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admits a family of steady states 

/'o\2\l/4 

M^) = ^7=^ ' ^ > 0' (5-25) 

all of which are unstable by Theorem 15.61 For the equation with rather than Krieger 
and Schlag [13] construct a stable manifold S associated with (fi. That is, if the initial data 

u{x, 0) = (pi{x) + i>o{x), dtu{x, 0) = i>i{x) 

correspond to a small perturbation {ipo,ilJi) G S, then the solution to f l5.24p exists globally and 
remains near the curve {v5a}a>o for all times. Numerical computations 0,(21] suggest that the 
stable manifold S ought to represent the borderline case between global existence and blow up 
of solutions. To see that Theorem 15.61 provides a partial proof of this conjecture, we note that 
the tangent plane to S at the origin is given by the equation 

a / xi^) i'oix) dx + / = 0, (5.26) 

where a > and — is the first eigenvalue of the linearized operator; see [13]. Thus the set of 
perturbations (12.21) to which our instability theorem applies is precisely the set of perturbations 
that lie strictly above the tangent plane to E at the origin. 
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